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£ ! Abstract 
i— » 

We consider an ensemble of large non-Hermitian random matrices of the form 
H + iA s , where H and A s are Hermitian statistically independent random 
N x N matrices. We demonstrate the existence of a new nontrivial regime 
of weak non-Hermiticity characterized by the condition that the average of 
[ — ' NTrA 2 is of the same order as that of TrH 2 when N — ► oo. We find explicitly 

the density of complex eigenvalues for this regime in the limit of infinite 
matrix dimension. The density determines the eigenvalue distribution in the 
crossover regime between random Hermitian matrices whose real eigenvalues 
are distributed according to the Wigner semi-circle law and random complex 
matrices whose eigenvalues are distributed in the complex plane according to 
S the so-called "elliptic law" . 

g 

o 

o 

Recently there was a growth of interest in statistics of complex eigenvalues of large 
^ ■ random matrices, both in physical and mathematical literature (see Refs. 1 - 20). 

Non-Hermitian random Hamiltonians appear naturally when one deals with quantum 
scattering problems in open chaotic systems |7|)|^, [lT|J15| ) [16|Jl8| , studies a motion of flux lines 
in superconductors with columnar defects [JHJ or is interested in chiral symmetry breaking 
in quantum chromodynamics ||20|| . Complex random matrices appear also in studies of dissi- 
pative quantum maps [^,^J. On the other hand, closely related ensembles of real asymmetric 
random matrices enjoy applications in neural network dynamics |5||12| and in the problem 
of the localization transition of random heteropolymer chains ||2 1|| . 



The actual progress in understanding the properties of random matrices with complex 
eigenvalues is rather limited, however. Most of the known results refer to the case of strong 
non-Hermiticity or asymmetry. Namely, they deal with those types of matrices for which 
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the real and imaginary parts of complex eigenvalues Zj = Xj + iYj are typically of the same 
order when the matrix dimension tends to infinity. 

As the most well-known example, we mention an ensemble of complex X x X matrices 
J with the density of joint probability distribution of matrix entries of the form 

V(J) = C N l exp [ - N Tr( J.P - r Re J 2 )} , (1) 

where r is real, -1 < r < 1, and C N = [n 2 Jo(l - r 2 )/X 2 ] N ^ 2 . In the large-iV limit the 
eigenvalues of J are uniformly distributed in the ellipse X 2 /(l + r) 2 + Y 2 /{I — r) 2 < Jq 
110 



Denote by Jki a matrix element of J. Eq. ([]]) implies that: (i) complex variables Jki are 
Gaussian with mean zero and variance (\Jki\ 2 ) = Jq/X ', (ii) the real and imaginary parts of 
J are statistically independent, and (iii) among all the random variables Jm, 1 < k, I < N, 
only Jm and are pairwise correlated with the magnitude of correlations determined by 
the covariance (JkiJik) — t Jq/N. The ensemble defined by Eq. (|l|) interpolates between the 
well-known Gaussian ensemble of Hermitian matrices (GUE) [F2^] and Ginibre's ensemble 
of complex matrices The degree of non-Hermiticity is controlled by the absolute value 
of r. If t = 1 or t = — 1 we have an ensemble of Hermitian or skew-Hermitian matrices 
and the above mentioned ellipse degenerates into an interval of real or purely imaginary 
eigenvalues. On the other hand, if r = then all Jki are mutually independent and we have 
a maximum of asymmetry. In that case the eigenvalues are uniformly distributed in the unit 
circle g],|22], 



The above mentioned ensemble of non-Hermitian random matrices can be represented 
in another form. Each matrix J can be decomposed into a sum of its Hermitian and skew- 
Hermitian parts: J = H + iA s , where H = (J + J t )/2 and A s = ( J - J t )/2z. By Eq. (fj), 
the Hermitian matrices H and A s are statistically independent and follow the Gaussian 
distributions V(H) oc exp{-2iV/[(l+r) J 2 } TrH 2 } andV(A s ) oc exp{-2JV/[(l-r) J 2 ] Tri 2 }, 
respectively. In this representation interpolation mentioned above becomes even clearer. For 
instance, when r approaches 1 variations of A around its mean value vanish and J = H in 
the limit of r = 1. 

For our purpose it is convenient to introduce new variables, J 2 = Jq(1 + r)/4 and 
v 2 = (1 — t)/(1 + r), and rewrite the ensemble defined by Eq. ([I]) as 

J = H + ivA, (2) 

where now each H and A is taken from identical and independent Gaussian ensembles of 
Hermitian matrices defined by the probability densities: 

V{H) oc exp ( - ^M 2 ) and V{A) oc exp ( - ^Tri 2 ) . (3) 

It is natural to use the parameter v 2 as a measure of non-Hermiticity (asymmetry) in 
our ensemble. As we mentioned, in the case of strong non-Hermiticity, i. e. when v 2 = 0(1) 
as X — > oo, the eigenvalues of J are uniformly distributed in an ellipse. 

In the present paper we demonstrate the existence of a new non-trivial regime of weak 
non-Hermiticity determined by the condition v 2 X = 0(1) as X —>■ oo. In other words, we 
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scale the parameter v with the matrix dimension N as v = a/y/N, a > 0. In that regime 
of weak non-Hermiticity the "elliptic law" is not valid any longer and should be replaced by 
another distribution. The main goal of our publication is to derive the explicit form of this 
distribution. 

In order to get better understanding of the chosen scaling it is instructive to consider 
first the limiting case of extremely weak non-Hermiticity. In this case one expects that the 
influence of adding the non-Hermitian matrices ivA to the Hermitian matrices H in Eq. 
can be treated by regular perturbation theory. 

When v — the eigenvalues Zj = X 3 of J = H are real and in the large-iV limit 
they are distributed according to the Wigner semi-circle law with the density p sc (X) = 
(27rJ) _1 y / 4 — (X/J) 2 . The mean separation A between adjacent eigenvalues can be esti- 
mated as A oc J/N. When a small non-Hermitian part ivA is taken into account, the 
eigenvalues Zj of J move into the complex plane and their shift can be estimated by first 
order perturbation theory. Namely, Re Zj = Xj and Im Zj = v((f)j \ A \<pj), where Xj and 
| 4>j) are the ordered eigenvalues of H and their normalized eigenvectors, respectively. Let 
us estimate the magnitude of this shift. Introducing the notation Pj for the projection on 
\<pj), one can write Im Zj = v TrAP,, so the variance of Yj = ImZj is 

where the angle brackets (-) H ^ denote the averaging over the random matrices H{A) ac- 
cording to Eq. (|3p. Since H is statistically independent of A, so is Pj. As a result, one can 
perform the averaging over A in Eq. (^) explicitly: 

4 T 2 

(( Y n A ) H =w {TTP ^^ (5) 

Now noticing that P 2 = Pj and TrP,- = 1 (as it must be for any projection), we obtain 

the deviation of Yj, ( Y 2 ) 1/2 = 2vJ/N. Since we are using perturbation theory, this de- 
viation has to be compared with the mean separation A of the unperturbed eigenvalues: 
( Y 2 ^^/A oc v\fN = a; so the case a < 1 corresponds to well-defined perturbation 
theory. 

The density p(X, Y) of complex eigenvalues Zj can be calculated easily in the perturba- 
tive regime a < 1. Indeed, p(X, Y) = (( N' 1 J2f=i S(X-X j )S(Y-Y j ) ) a )h- Now employing 
the Fourier representation for S(Y), one finds that 

= J- £ r dke ™( s(x - Xj ) exp [ - WM-am ) 

2irN ^J-oo \ v L 2 \ /h 

4EU(X-^)4X>,e-exp[-i(^) 2 ] 
= p sc (X)p(Y), 
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where 



P(Y) 



2-na 



- exp 



y2 
2d 2 



(6) 



and a = a J/N. We see that for a < 1 the density of eigenvalues in the complex plane X+iY 
is a simple product of the Wigner semicircular density p sc (X) = (2ix J)~ l ^J<i — (X/ J) 2 and 
the Gaussian p(Y) from Eq. (^). On the other hand, it is natural to expect that when a ^> 1 
p(X, Y) goes back to the mentioned uniform distribution inside an ellipse. 

In order to get access to the distribution of eigenvalues in the complex plane E = X + iY 
in the case of a = 0(1) we use the fact that the two-dimensional density of these eigenvalues 
p(X, Y) can be found if one knows the "potential" |§: 



2tcN 



In Det[(E - J)(E - J) t + k 2 I] 



A/H 



(7) 



in view of the relation: p(X,Y) = lim K ^o d 2 <&(X, Y, k), where d 2 stands for the two- 
dimensional Laplacian d\ + d\. To determine the potential for the present case of almost- 
Hermitian Gaussian random matrices we follow the method suggested earlier by two of us 
|T8| and restore the potential from its derivative: 



<9 2 $ 
dn 2 

Z(K b , K) 



d 



lim 







((In Z(K b ,K)) A ) h 



2tcN dK "t^K 8k " s ~ v ' "' /A/H 
Det[(E - J)(E - Jy + k 2 I] 
Det[(E - J)(E - J)t + k 2 J]' 



(8) 



It is convenient for our purpose to write down the determinants in the denominator and 
numerator of the generating function as: 



Det[(E -J)(E- J) f + k 2 J] = Det 



M h 



kJ i(E — J) 
i{E - jy kJ 



M b l 



(9) 



Det[(E -J)(E- Jy + k 2 I] = Det 



where we used two matrices: M b 



M 



i 

v/2 



i(e - Jy kI 

kI i(E - J) 



M' 1 



(10) 



/ -il 

i ii 



M, 



M b l and M 



i 

V2 



I -I 

i i 



Then one can find the following convergent representation for the generating function 
Z(ku, k) in terms of a Gaussian integral over both commuting and anticommuting (Grass- 
mannian) variables: 



[-l) N Z(K b , K ) = J [d$] exp{-£ (*) - 



c (y) = K b (&kLV) + iX(m ] Lm) - Y(&& y) + (k- K b )(m ] km) 



;n) 



d(m) = -itf(H®L)V 



a 



#1/2 
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where ^ 



[dM] = U. p=lt2 dSpdS^dx P dxl, with S p and x P being 
iV— component vectors of complex commuting and Grassmannian variables, respectively. 
The 4x4 matrices A, L, ctq and K are (block) diagonal of the following structure: 



A 
K 



diag(l, -1,1,-1); L = diag(l, -1, 1, 1); 



diag(0, 0,1,-1); cf 



diag(zS :c , SJ; 




and T, x = 

Due to the normalization condition Z(n, k) = 1 it is enough for our purpose to calculate 
the average ({Z(k,Ki)) a ) h . This can be done easily for the Gaussian distributions Eq. @ 
and leads to the terms quartic in the components of the supervector \1/ in the exponent of the 



integrand, Eq. (|TT|) . Further evaluation goes along lines suggested by Efetov in the theory 
of disordered systems P5| . A detailed introduction to the method as applied to random 
Hermitian matrices can be found in the review [f24 |. Here we outline only the general 
strategy: 1) To decouple quartic terms in the exponent of the integrand of the generating 
function by introducing a set of auxiliary integrations (the so-called Hubbard- Stratonovich 
transformation, see |24|); 2) To integrate out the ^/-variables explicitly; and 3) Exploiting 
the limit N — ► oo to integrate out some ("massive") degrees of freedom in the saddle-point 
approximation. After this is done the integral over the remaining ("massless") degrees of 
freedom can be represented in a form of the so-called zero- dimensional graded nonlinear 
u — model introduced into physics by Efetov [23 . 

After this set of standard manipulations one arrives at the following expression for the 
density px{v) of the scaled imaginary parts yj = 2irp sc {X)YjN of those eigenvalues Zj = 
Xj+iYj whose real parts Xj fall within a narrow window around the point X of the spectrum: 



1 d 2 f°° 
Px(y) = 2dfJo 



duu(f>(y, u) 



(12) 



where 



C) r A A r 111 A A 111 A A A "I 

</>(?/, u) = i ir dfi(Q)Sti(KQ) x exp - — Str(QA) - ^Str(Qa) - Str(QaQa) , (13) 
ou J 1 2 2 16 J 

where a = 27rp sc (X)a, a = diag(S a; , T, x ) and the (graded) matrices Q satisfying Q 2 = — 1 
are taken from the graded coset space {7(1, l/2)/{7(l/l) £g> {7(1/, 1), whose explicit parame- 
terization can be found in [23,24|. We also use the notation Str for the graded trace. 



Still, the evaluation of the integral over the graded coset space in Eq. flT3p and subsequent 
restoration of the density px(y) is quite an elaborate task. We skip intermediate steps 
and present the final expression for <p(y, u) which one obtains after integrating out the 
Grassmannian variables: 



I(zi, z 2 ) 
where 



(I(zi,z 2 )); 



1 roo roo 



2ll j — oo j —oo 

a 2 



dz\dz 2 exp 



'A + z\) 



Hzi,z 2 ) 



(14) 



y [F 2 d y F x - F 1 d v F 2 \ - - F,F 2 + d y F l d y F 2 + d u F 1 d u F 2 - F 2 d l y F x - F x d l y F 2 
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e -\J u2 +vl sinh Ju 2 + yj 
Fi(y, u, zi, z 2 ) = . ; F 2 (y, u, z x , z 2 ) = . v ; (15) 

V^ 2 + vi v n + $ 

and yb = y ~ azi/ \[2\ y f = y- iaz 2 / y/2. 

It is easy to check that both functions i*i and F 2 satisfy the relation: 

dlF^ = F h2 - -d u (ud u F 1>2 ) . (16) 

Now substitute Eq. ([L4]) into Eq. ( fl2|) and notice (by multiple exploitation of Eq. (|TjJ) ) 
that I(zx, z 2 ) can be written as a full derivative over u. As a result, one obtains: 



px{y) = {[2 + yd y ] [u{F l d u F 2 - F 2 d u F x )\ | u ^ ) z - (17) 

where {■■■)z denotes the integration over z x ,z 2 with the Gaussian weight, see Eq. (|14|). 
When performing the limiting procedure u — > in Eq. (|T7|). one notices that all terms are 
proportional to a 5-function: 5(yb), or its derivatives. This fact allows one to perform the 
integration over z\ explicitly, and after some simple algebraic manipulations one arrives at 
the following expression: 

p x {y) = 77=^exp f-^J ^ 1 ^cosh(^)exp(-a 2 t 2 /2), (18) 



where a = 2irp sc (X)a and y = 2irp sc (X)YN, X and Y being the real and imaginary parts 
of the complex eigenvalues in the ensemble of random matrices J given by Eqs. (@)-(0)- 

The distribution Eq. (|18|) constitutes the main result of the present publication. It 
correctly reproduces all the anticipated limiting cases. When a ^> 1 one can effectively 
put the upper boundary of integration in Eq. (|18D to be infinity due to the Gaussian cut- 
off of the integrand at f ~ 1/a C 1. This immediately results in the uniform density 
Px{y) — (2a 2 )" 1 inside the interval \y\ < a 2 and vanishing density outside this interval. 
Recalling the definition of the variable y and the parameter a, we are back to the familiar 
"elliptic law" . In the opposite limiting case a<l one obtains: 

1 / y 2 n sinh y 
Px{y) = , u/2 exp ( - — ) 

which matches the perturbative result Eq. (^) as long as y 1. The fact that remote 
tails y ^ 1 are not captured correctly by perturbative treatment can be understood easily: 
the condition y ~ 1 means that the imaginary parts Yj of the complex eigenvalues is of 
the same order as the mean spacing A = (p sc (X)N)~ l between the real eigenvalues of the 
unperturbed Hermitian matrix H. First order perturbation theory is obviously insufficient 
to obtain those complex eigenvalues correctly. 

It is worth mentioning that we have restricted ourselves to Gaussian distributions of 
random matrices H and A mainly for the sake of simplicity of presentation of the supersym- 
metry calculations. More general distributions of matrices H and A can be considered and 
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FIG. 1. Histograms of distributions of the scaled imaginary parts of complex eigenvalues for a) 
Gaussian and b) Bernoulli distribution of matrix elements, see the text. 



super symmetry calculations can be performed along the line of [^,^|. For instance one 
can consider random matrices H and A whose entries are independent Bernoulli variables 
taking ±6 with equal probability or are uniformly distributed in an interval. Actually, we 
have arguments that the universality of the eigenvalue distribution in the crossover regime 
is extremely high and is not restricted to the case when the corresponding Hermitian matri- 
ces H have the semi-circle eigenvalue density. For example, one can consider non-Gaussian 
distributions of matrices H and A with the density of the form V(H) oc exp{— NTtV(H)} 
(the so-called invariant ensembles, see f27H33|). Then the formula Eq. (|T8| ) is still valid, 
provided the semi-circular density p sc in the definitions of y and a is replaced by the actual 
mean level density at given point of the spectrumF] . 

The distribution Eq. flT8| ) has been derived in the limit of infinite matrix dimension. 
However, numerical computations show that Eq. ([18]) well approximates the distribution of 
complex eigenvalues for ensembles of random matrices (§) of a moderate dimension. In Fig. 1 
we present results of numerical diagonalization of random matrices of dimension N = 30. 
For each plot we generated a set of 20000 random matrices J = H + i^jjjA sampling a 
Gaussian distribution (random matrices follow distributions of Eq. @ with 2J 2 = 1) for 
plot (a) and a Bernoulli distribution (functionally independent matrix elements of H and 
A are statistically independent random variables taking ±l/\/2iV with equal probability) 
for plot (b). Then matrices were diagonalized and their eigenvalues Zj falling into a small 
energy window, |Re Zj\ < 0.2, were selected. The imaginary parts of the selected eigenvalues 
were counted and the corresponding histograms along with the infinite- N theoretical density 
Px=o(y) (solid curves) were plotted against the scaled variable yj = 2irp sc (0)YjN. In both 
plots, the ordinate is scaled in units of max y px=o(y)- 



1 We assume that the matrices H, A normalized in a way ensuring (TrH 2 ) = (Tr^ 2 ) = J 2 N. 
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It is worth mentioning that our derivation based on the supersymmetry formalism p3|j24 



being in general quite satisfactory should be considered as a heuristic one from the point of 
view of rigorous mathematics. In a more detailed publication |34] we shall give a rigorous 
mathematical derivation of our main result, Eq. fll8|) for the case of Gaussian Matrices H, A 
and also show the mentioned universality of the found distribution for several classes of more 
general random matrices. 

At the moment we failed to derive the analogous expression for almost-symmetric 
(slightly asymmetric) real random matrices due to unsurmountable technical problems. At 
the same time there are reasons to suspect that for the latter case the distribution might be 
different in its form from Eq. flT8|). It would be of much interest to try to attack this problem 
by different methods, e.g. starting from the known joint probability density for eigenvalues 
of asymmetric random matrices |iO| , |T3|| . 

In conclusion, we derived the explicit expression describing the density of complex eigen- 
values of almost-Hermitian random matrices. It describes the crossover regime from the 
Wigner semi-circle law typical for random Hermitian matrices to the "elliptic law" typical 
for strongly non-Hermitian matrices. 

The authors acknowledge financial support of Deutsche Forschungsgemeinschaft under 
Grant No SFB237. Numerical computations were performed at Institut fiir Mathematik, 
Ruhr-Universitat-Bochum with the help of the MATLAB software package. 
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